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ABSTRACT 

Genetic  variation  in  the  plasma  enzyme  butyryl cholinesterase  (BCHE)  affects  the  response  of 
humans  to  xenobiotic  agents.  Mutations  in  BCHE  are  responsible  for  the  majority  of  cases  of 
prolonged  apnea  following  the  administration  of  sued nyl choline.  In  addition,  genetic  variation 
in  BCHE  is  linked  to  sensitivity  to  both  organophosphate  and  carbamate  compounds,  including 
the  deployment  related  drugs  pyridostigmine,  physostigmine,  heptyl  physostigmine  and  SDZ- 
ENZ  713.  The  study  described  in  this  report  successfully  screened  the  four  coding  regions  (and 
surrounding  intronic  regions)  of  BCHE  for  both  novel  and  known  polymorphisms.  High 
Resolution  Melt  (HRM)  and  sequencing  procedures  revealed  12  different  genetic 
polymorphisms,  and  35/  51  individuals  were  shown  to  carry  at  least  one  BCHE  polymorphism. 
Eight  of  the  polymorphisms  had  been  documented  previously,  two  of  which  have  been 
reported  in  individuals  with  sued  nyl  choline  and  pyridostigmine  sensitivities.  Of  the  four  novel 
polymorphisms  found,  two  are  predicted  to  change  an  amino  acid  in  BCHE.  This  study  has 
demonstrated  clearly  that  unique  and  important  functional  genetic  variation  may  be  found  and 
there  may  be  a  small  subset  of  individuals  in  the  ADF  with  enhanced  sensitivity  to 
sued  nyl  choline  and  and  organophosphate  and  carbamate  compounds. 
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The  Identification  of  Butyryl cholinesterase  (BCH  E) 
Polymorphisms  in  a  Small  Australian  Defence  Force 

Cohort 


Executi  ve  Su  mmary 


Deployed  military  personnel  may  be  exposed  to  a  wide  range  of  chemicals  and  drugs. 
Following  the  Gulf  War,  some  veterans  reported  illnesses  that  might  have  been  the 
result  of  chemical  exposures,  and  significant  correlations  between  perceived  illness  and 
chemical  exposure  have  been  reported.  Reported  chemical  exposures  include  the  insect 
repellent  N,  N-diethyl-m-toluamide  (DEET),  insecticides  such  as  permethrin  and 
chlorpyrifos  used  to  protect  against  insect  borne  diseases  and  the  neuro-protective 
agent  pyridostigmine  bromide  used  to  protect  against  nerve  agent  attack. 

There  is  a  wide  range  of  chemicals  and  drugs  that  interact  with  the  cholinergic  nervous 
system  including:  nerve  agents,  nerve  agent  p re- treatments,  insecticides  and 
therapeutic  drugs.  Sensitivity  to  cholinesterase  inhibitors  and  effective  nerve  agent  pre¬ 
treatment  rely  on  the  activity  of  the  enzymes  acetylcholinesterase  (ACHE), 
butyryl  cholinesterase  (BCH  E)  and  paraoxonase  1  (PON  1),  and  these  genes  are  subject 
to  genetic  variation.  Genetic  variations  in  these  genes  have  been  shown  to  alter  the 
level  and  activity  of  their  enzyme  products,  and  effectivedrug  treatment  may  therefore 
be  i  nfl  uenced  by  the  genotype  of  the  i  ndi  vi  dual . 

BCHE  hydrolyses  chemicals  containing  ester  bonds  such  as:  drugs  acting  at  the 
neuromuscular  junction  (sued nyl choline),  local  anaesthetics  (procain,  chloroprocaine 
and  cocaine)  and  heroin.  It  is  also  a  biological  scavenger  against  organophosphorus 
and  carbamate  compounds  used  as  pesticides  and  nerve  agents.  Individuals  carrying 
heterozygous  or  homozygous  inactivating  alleles  have  been  shown  to  have  a  lower 
capacity  to  interact  with  and  detoxify  drugs.  Variation  in  an  individual's  genetics  may 
therefore  result  in  failure  to  achieve  a  desired  therapeutic  effect  with  nerve  agent  pre- 
and  post-treatments  or  development  of  adverse  drug  reactions.  In  a  military  setting, 
genetic  variation  in  BCH  E  may  result  in  altered  response  to  deployment  related  drugs 
including  pyridostigmine,  physostigmine,  heptyl  physostigmineand  SDZ-ENZ  713. 

Our  initial  work  has  focused  on  characterisation  of  the  genetic  variation  in  BCHE  in  a 
small  cohort  (n=51)  of  ADF  personnel.  Using  the  LightCyder  480  system,  a  protocol 
was  established  using  High-Resolution  Melting  Curve  Analysis  to  screen  the  BCHE 
gene.  Any  PCR  amplicons  displaying  altered  melting  profiles  were  subjected  to  direct 
sequencing,  and  variants  were  identified. 
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This  study  has  identified,  in  a  subset  of  ADF  members,  genetic  variation  in  BCHE 
which  may  confer  increased  sensitivity  to  a  range  of  chemicals  and  drugs  of  interest  to 
Defence.  Initial  screening  has  identified  12  different  polymorphisms,  8  known  and  4 
novel.  Two  of  the  known  polymorphisms  have  known  functional  effects  with 
individuals  carrying  these  polymorphisms  exhibiting  abnormal  sensitivity  to 
sued  nyl  choline. 

Future  work  will  involve  the  characterisation  of  the  four  novel  variants  to  determineif 
they  are  associated  with  a  functional  effect  on  either  protein  levels  or  enzyme  activity. 

In  summary,  this  work  demonstrates  clearly  the  existence  in  the  ADF  of  genetic 
variations  that  may  have  a  significant  impact  on  the  response  of  individuals  to 
pharmaceuticals  and  their  sensitivity  to  chemical  warfare  agents.  Furthermore,  these 
susceptible  individuals  may  be  identified  using  a  simple,  non-invasivetest  procedure. 
It  is  also  likely  that  functional  variants  unique  to  Australian  populations  exist  and  that 
these  will  need  to  be  identified  in  order  to  accurately  predict  chemical  and  drug 
response. 
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1.  Introduction 

Organophosphorus  compounds  (O  Ps)  are tri  esters  of  phosphori  c  aci  d  and  thei  r  major  use  i  s  as 
insecticides.  Chemical  warfare  agents  that  affect  the  nervous  system  and  nerve  agents  (eg. 
soman  and  sarin)  also  belong  to  this  class  of  chemicals.  On  acute  exposure,  OPs  cause 
neurotoxicity  in  the  form  of  a  cholinergic  syndrome  (i.e.  overstimulation  of  muscarinic  and 
nicotinic  acetyl  choline  receptors  in  the  central  and  peripheral  nervous  systems)  resulting  from 
accumulation  of  acetylcholine  in  the  synaptic  cleft  as  a  consequence  of  inhibition  of 
acetylcholinesterase  (AChE).  In  addition,  both  acute  and  repeated  low-level  exposure  to  some 
OPs  is  associated  with  delayed  neuropathy  and  can  result  in  long  term  cognitive  effects  [1, 2], 

Human  butyryl  cholinesterase  (6  CHE)  is  an  enzyme  found  in  plasma  and  many  other  parts  of 
the  body.  The  enzyme  hydrolyses  drugs  containing  ester  bonds  such  as  drugs  acting  at  the 
neuromuscular  junction  (sued nyl choline),  local  anaesthetics  (procaine,  chloroprocaine  and 
cocaine)  and  heroin  [3].  It  is  also  a  biological  scavenger  for  organophosphorus  and  carbamate 
compounds  used  as  pesticides  [4].  The  rapid  hydrolysisof  succinylcholine  by  BCWE  accounts 
for  its  short  duration  of  action.  Mutations  in  the  BCHE  gene  affecting  the  function  of  this 
enzyme  can  I  ead  to  prol  onged  neuromuscul  ar  bl  ockade  by  sued  nyl  chol  i  ne.  I  n  a  mi  I  itary  setti  ng, 
abnormal  BCHE  enzyme  levels  may  result  in  altered  response  to  deployment  related  drugs 
including  pyridostigmine,  physostigmine,  heptyl  physostigmineand  SDZ-ENZ  713  [5, 6], 

The  B  CHE  gene  is  located  on  chromosome  3q26.1-q.  26. 2.  The  genomic  region  for  BC  HE  spans 
over  70  kb  and  has  4  exons  and  3largeintrons.  The  gene  encodes  a  transcript  of  2400bp  in  size 
and  a  protein  of  574  ami  no  acids.  To  date,  more  than  100  polymorphisms  have  been  identified; 
however  few  havebeen  studied  fully.  In  general,  BCHE  polymorphisms  have  been  shown  to 
produce  enzymes  w  ith  varyi  ng  I  evel  s  of  catal  yti  c  activity.  The  mol  ecul  ar  bases  of  several  geneti  c 
variants  of  BC  HE  havebeen  reported,  such  astheAtypical  gene, fluoride-resistant  gene,  silent 
gene,  K  variant,  J  variant  and  H  variant  [7-10].  In  addition  there  are  a  number  of  additional 
BCHE  polymorphisms  which  result  in  a  protein  with  no  enzymatic  activity.  Enzymes  with 
activity  below  10%  of  the  wild  type  enzyme  are  cal  led  "silent"  variants  [11],  An  example  of  the 
rangeof  polymorphisms  identified  within  BCHE  are  summarised  in  Appendix  A. 

Currently,  littlegenetic  analysis  of  BCHE  has  been  performed  in  samples  from  an  Australian 
population.  A  single  study  performed  in  2003  screened  the  BCHE  gene  of  65  individuals 
suffering  prolonged  post-succi  nyl  choline  apnea.  Overall  the  Atypical  and  K  variant  mutations 
were  highly  prevalent,  however,  six  new  BCHE  mutations  were  identified,  and  itwas  found 
most  individuals  carried  multiple  homozygous  mutations  [12].  This  demonstrated  unique 
heterogeneity  in  a  proportion  of  the  Australian  population  and  that  novel  functional  variants 
exi  st  i  n  th i  s  popu  I  ati  on .  There  i  s  I  i  ttl  e  i  nformati  on  avai  I  abl  e  on  the  I  evel  of  geneti  c  vari  ati  on  i  n 
theAustralian  population,  and  the  I  evel  of  inter-individual  variation  i  n  the  expressi  on  of  BC  HE 
is  unclear. 

Therefore,  combinations  of  phenotypes  and  expression  levels  of  BCHE  involved  in  OP 
detoxification  might  be  expected  to  occur  in  military  personnel  and  define  a  subgroup  who  are 
potentially  more  susceptible  to  adversedrug  effects  than  the  majority  of  ADF  members. 

Our  approach  to  BCHE  genotyping  considered  issues  of  time-  and  cost  effectiveness. 
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The  ai  m  of  thi  s  project  was  to; 

•  Screen  for  novel  sequence  variations  within  thecoding  regions  of  the  BCHE  gene. 

•  I  dentify  the  presenceof  any  common  polymorphisms  of  functional  relevancetotheADF. 

•  Comparethegenotypefrequency  of  these  polymorphisms  and  compareto  the  level  and 
types  of  geneti  c  vari  ati  on  w  i  th  that  found  i  n  other  popul  ati  ons  (U  S,  Eu  ropean,  Chi  nese). 


2.  Material  and  Methods 


2.1  Population 

The  experimental  procedureswereapproved  by  the  Australian  DefenceH  uman  Research  Ethics 
Committee(ADHREC  protocol  418 /  05),  and  all  guidelines  were  adhered  to.  Participants  were 
51  Australian  Defence  Force  soldiers  who  provided  10  ml  of  whole  blood  as  part  of  a  larger 
medical  screening  process.  All  samples  were  unlabel  led  at  the  time  of  col  lection  to  allow  for  a 
double-blind  analysis,  and  no  participant  could  be  identified  from  their  biological  material. 
Samples  used  as  controls  were  obtained  from  the  Red  Cross,  agreement  number  08- 11V I C-05. 

2.2  Preparation  of  G  enomic  D  N  A 

Genomic  DNA  was  isolated  from  whole  blood  using  the  QIAamp®  DNA  Blood  Mini  kit 
(Qiagen),  using  "Blood  or  Body  Fluid  Spin  Protocol"  withthefollowing  modifications.  A  total  of 
200 pi  of  whole  blood  was  used  in  each  extraction  and  samples  were  eluted  in  200 pi  of 
RNAase/  DNAse  free  water  and  stored  at4°C. 

2.3  High  Resolution  Melt  Strategy 

All  primers  were  designed  using  Primer3  software  which  is  freely  available  at 
http:/  /  frodo.wi.mit.edu/  primer 3/  input.htm.  The  following  parameters  were  applied  to  all 
primers;  primer  length  18-27  bp,  GC%  =  30-70%,  Tm  =  59-60°C  and  amplicon  size  =  100- 
250  bp.  Amplicons  were  designed  to  cover  all  exonic  regions  of  the  BCHE  gene  and  flanking 
intronic  regions.  Primer  sets  were  designed  so  that  multi  pie  amplicons  across  an  exon  had  at 
least  a  30  bp  overlap.  Primer  sequences  and  amplicon  size  are  described  in  Technical  report 
DSTO-TR-2498.  Figure  1  describes  the  position  of  the  24  amplicons  within  the  BCHE  gene 
structure. 
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6272  bp 


exon  1  exon  2 


159  bp  1525  bp 


exon  3  exon  4 

167  bp  603  bp 


Exonl  :1  amplicon 
Exon  2  :13  amplieons 
Exon  3  :4  amplieons 
Exon  4 :5  amplieons 
3'  UTR:  1  amplicon 


Figurel:  Schematic  representation  of  the  gene  structure  and  relative  position  of  amplieons  used  in 
analysis  BCHE.  BC HE  consists  of  4  exons,  exon  1  is  159  bp,  exon  2  is  1525  bp,  exon  3  is 
161  bp  and  exon  4  is  603bp.  BCHE  messenger  RN A  transcript  is  2454  bp  is  size  and  the 
initiating  AT  G  is  I  ocated  at  base  9  of  exon  2  (red)  and  stop  codon  is  at  posi  ti  on  1 25  bp  of  exon 
4.  In  total,  exon  1  was  covered  by  1  amplicon  (blue),  exon  2;  13  amplieons  (green),  exon  3;  4 
amplieons  (aqua),  exon  4;  5  amplieons  (pink)  and  3'UTR  1  amplicon,  purple. 


Real  -ti  me  PC  R  ampl  ifi  cati  on  of  theB  C H  E  genefragment  was  performed  i  n  a  L  i  ghtCyd  er®  480 
H  i  gh  Resol  uti  on  M  el  ti  ng  M  aster  mi  x  contai  ni  ng  a  saturati  ng  H  RM  dye  (Roche  D  i  agnosti  cs  Co, 
Mannheim,  Germany).  Genomic  DNA  from  unknown  samples  (5ng)  was  added  to  a  20 pi 
reaction  mixturecontaining0.2  pM  of  each  primer  and  lxHigh  Resolution  Master  mix.  MgC  1 2 
concentrations  varied  according  to  the  primer  pair  and  are  described  in  technical  report 
DSTO-TR-2498.  Samples  were  also  spiked  with  a  20%  homozygous  control  sampleas  described 
in  DSTO-TR-2498.  PCR  was  carried  out  using  a  touchdown  protocol,  with  annealing 
temperatures  ranging  from  62°C  to  53°C.  H  igh-resolution  melting  curvedatawereobtained  at  a 
rate  of  25  acquisitions  per°C.  PCR  conditions  are  described  in  Table! 


Tablel:  LightCycler  480  Reaction  conditions 


Target  (°C) 

Acquisition 

Mode 

Hold 

(hh:mm:ss) 

Ramp  Rate 
(°C/s) 

Acquisitions 

(per°C) 

Pre-incubation 

95 

none 

00:10:00 

4.4 

0 

Amplification  45 cycles 

95 

none 

00:00:10 

4.4 

0 

62  -  53* 

none 

00:00:10 

2.2 

0 

72 

single 

00:00:10 

4.4 

0 

High  Resolution  Melting 

95 

none 

00:01:00 

4.4 

0 

40 

none 

00:01:00 

2.2 

0 

65 

none 

00:00:01 

4.4 

0 

95 

continuous 

0.02 

25 

cooling 

40 

none 

00:00:30 

2.2 

0 

*  53°C  is  second  temperature  target,  with  a  step  size  of  0.5°C  per  cycle 
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Melting  curves  were  analysed  using  the  LightCyder  480  software.  Briefly,  melting  curves  for 
individual  amplification  productswerenormalised  using  a  pre-melt  signal  setto  a  relative  value 
of  100%  and  post  melt  signals  set  to  a  relative  value  of  0%.  The  melting  shift  was  calculated 
automatically  by  the  software  using  a  default  adjustment  value  of  5%  for  all  analyses.  Usinga 
cu  rve  matchi  ng  al  gorithm,  sampl  es  were  cl  ustered  i  nto  groups  and  difference  pi  ots  exami  ned . 
Thesensitivity  level  was  set  on  0.3  or  0.6for  a  high  sensitivity  setting.  Higher  sensitivity  settings 
identify  smaller  fluorescent  differences  which  are  required  when  multiple  melting  domains 
exist  within  theamplicon  examined.  Amplicons  displaying  alternative  melting  profiles  were 
confirmed  by  sequencing  an  independent  PCR  product. 

2.4  DNA  sequencing 

PCR  amplifications  were  performed  in  a  20  pi  reaction  containing  0.5  pM  of  each  primer,  5  ng 
gDNA,0.8U  FastStartTaq  (Roche),  0.5  mM  of  each  dNTP  (Roche),  2.0mM  MgCband  lx  PCR 
buffer  (Roche),  using  amplification  conditions  described  in  Table  3.  Samples  were  analysed  on 
1.2%  Agarose  gels  to  confirm  single  products  of  correct  size  and  subsequently  purified  using 
EXO-SAPIT  (GE  Healthcare).  Samples  were  sequenced  using  the  BigDye  Terminator  kitV3.1 
(Applied  Biosystems)  and  cleaned  using  Xterminator  (Applied  Biosystems)  according  to  the 
manufacturer's  instructions.  Analysis  of  the  sequencing  reactions  was  performed  on  an 
A  Bl  3130x1  (A  ppl  i  ed  Bi  osystems).  Sequence tracefi  I  es  were  anal  ysed  for  SN  Ps  usi  ng  Sequencher 
(GeneCodes,  Genesearch). 


3.  Results 

In  this  study,  we  tested  the  ability  of  HRM  using  the  LightCyder  480  to  screen  for  DNA 
vari  ati  ons  i  n  the  cod  i  ng  regi  ons  of  theB  C  H  E  gene.  I  n  total ,  24 fragments  i  n  each  of  the  51  gD  N  A 
sampl  es  were  scanned . 

Figure  2  shows  the  melting  curves  and  subtractive  difference  plots  for  all  exons  of  the BCHE 
gene  i  n  w  hi  ch  at  I  east  one  vari  ati  on  w  as  i  d  entifi  ed .  The  mel  ti  ng  cu  rves  i  d  enti  f  i  ed  D  N  A  sampl  es 
with  possible  polymorphisms.  Subtractive  difference  plots  were  obtained  by  subtracting  the 
curve  of  the  sample  of  interest  from  that  of  a  control  sample  with  normal  DNA  sequence. 
Abnormalities  initially  identified  in  the  melting  curves  were  confirmed  and  magnified  by  the 
subtractive  difference  plot.  These  samples  were  then  amplified  by  PCR  and  sequenced 
(Appendix  B). 
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Exon  1  fragment  1 


Exon  2  fragment  3 


Exon  2  fragment  4 


^  Novel  3 

(C/T) 

Novel  3 

1 

1  /  — 

4 

Exon  2  fragment  6 


Exon  2  fragment  7.2 


rsl63497C0  (C/G) 

ii 

I. 

R  efer  to  figure  legend  page  6 
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Exon  3/  Intron  2  fragment  11 


rs!803274(G/A) 


Exon  4  fragment  14 


3'UTR  fragment  15 


Figure2:  M  elting  (left)  and  subtractive  difference  plots  (right)  of  normal  and  abnormal  fragments  of  the 
BCHE  gene.  Note:  theshifting of theplotsisduetotheindicated mutations.  Fragmentswere 
analysed  by  PCR  amplification  with  theHRM  dye  followed  by  a  high-resolution  thermal 
den  atu  rati  on.  Subtractive  plots  on  the  right  were  obtained  by  subtracting  each  curve  from  a 
non-polymorphiccontrol.  N  ote:  each  polymorphism  produced  a  uniquemelting  and  difference 
plot 


DNA  variations  in  the  B  CHE  gene  identified  in  the  51  subjects  are  summarised  inTable2. 

Of  the  51  subjects,  35  had  at  I  east  one  SN  P  i  n  theS  CH  f  regi  ons  exami  ned  .Two  polymorphi  sms 
are  well  described,  K  variant  and  Atypical  variant.  These  alleles  both  have  altered  function 
compared  to  wild-type.  K  variant  SNP  results  in  a  protein  with  30%  reduced  activity  and 
Atypical  is  a  silent  allele.  In  the  case  of  K  variant,  2  samples;  SOTG8-68and  SOTG8-76were 
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homozygous  for  the  mutant  allele.  Samples  RAN  1-12  and  RAN  1-21  were  shown  to  be 
heterozygous  for  both  polymorphisms.  Sequencing  analysis  uncovered  4  new  BCHE  mutations, 
not  currently  in  public  databases.  Two  of  theselle384Val  and  Ala212Val  result  in  amino  acid 
substitutions  in  exon  2.  In  summary  12  different  polymorphisms  were  identified  and  these  are 
su  mmar i  sed  i  n  Tabl  e  2. 


Table2:  Summary  of  DNA  variation  in  BCHE  gene,  (gDNA  NCBI  Reference  Sequence: 

NG_009031.1;  cDNA  NCBI  R eference Sequence:  N M  _000055. 2,  protein  N P_000046.1m 
M  ap  to  genome  build  37.1)  identified  in  HRM  analysisof51  subjects. 


Position  in  gene 

Amplicon 

Name  of  SNP 

DNA  variation 

Number  of  positive 
samples 

Exon  1 

5'  UTR 

1 

rsll26680 

g.5129G>A 

c.-32G>A 

G/  G  45/  51 

G/  A  6/  51 

exon  2 

3 

rsl799807 
Atypical  variant 

g,11725A>G 
C.293A  >G 
p.Asp98Gly 

A/  A  49/  51 

A/  G  2/  51 

4  and  5 

N  ovel  3 

g. 12067  C>T 
c.  801  C>T 
p.  Ala212Val 

C/  C  50/  51 

C/  T  1/  51 

6 

rsl6849700 

g,12281G>C 

c.849G>C 

p.Glu283Asp 

G/  G  49/  51 

G/  C  2/  51 

7.1  and  7.2 

Novel  1 

g. 15509  C>T 
c,1244C>T 
p.  Val384Val 

C/  C  49/  51 

C/  T  1/  51 

7.2  and  8.1 

N  ovel  2 

g.l2582A>G 
C.1317A>G 
p.l  le384Val 

A/  A  50/  51 

A/  G  V  51 

exon  3  (i  ntron  2) 

11 

rs55781031 

g.56034T>C 

c.1518-121T>C 

T/  T  46/  51 

T/  C  5/  51 

exon  4 

13, 13.1, 13.2 

rsl803274 

K  variant 

g.68974G>A 

C.1699G>A 

p.Ala567Thr 

G/  G  29/  51 

G/  A  21/  51 

A/  A  2/  51 

exon  4 

3'UTR 

14 

rs3495 

g.69273G>A 

c.2159G>A 

G/  G  27/  51 

G/  A  24/  51 

3'UTR 

15 

rs55849903 

g.69433C>T 

c.2225C>T 

p.N/A 

C/  C  49/  51 

C/  T  2/  51 

3'UTR 

15 

Novel  4 

g.C>T  69474 
c.2366C>A 
p.N/A 

C/  C  50/  51 

C/  T  1/  51 

near  gene  3' 

15 

rs77307863 

g.71985803C>A 

C/  C  42/  51 

C/  A  8/  51 

A/  A  V  51 

Of  the  samples  with  polymorphisms,  (24/  35)  had  compound  gene  defects  (i.e.  compound 
heterozygotes  and  homozygotes)  and  this  is  summarised  in  Table  3.  Five  samples  RAN  1-18, 
RANI-9,  RANI-4,  RANl-land  SOTG  8-45  identified  4  SNPs  within  the  regions  examined. 
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Table3:  Surrmary  of  polymorphisms  found  in  BCHE  in  the  regions  examined  fcr  each  individual 


Sample 

Exon  1 

Exon  2 

intron  2 (near  3'  ) 

Exon  3 

Exon  4 

3UTR 

SOTG8-45 

NP 

rsl6849700(G/C) 

NP 

NP 

rsl803274  (A/  G),rs3495  (G/  A) 

rs55849903(C/T) 

SOTG8-46 

NP 

NP 

NP 

NP 

NP 

novel  4G/  A 

SOTG8-47 

NP 

NP 

NP 

NP 

rsl803274 (A/  G),  rs3495  (G/  A) 

NP 

SOTG8-48 

NP 

NP 

NP 

NP 

NP 

rs77307863  (G/  T) 

SOTG8-49 

NP 

NP 

NP 

NP 

rsl803274 (A/  G),  rs3495  (G/  A) 

SOTG8-50 

NP 

NP 

NP 

NP 

rsl803274 (A/  G),  rs3495  (G/  A) 

NP 

SOTG8-51 

NP 

NP 

NP 

NP 

NP 

NP 

SOTG8-52 

NP 

NP 

NP 

NP 

NP 

rs77307863  (T/  T) 

SOTG8-53 

NP 

NP 

NP 

NP 

rsl803274 (A/  G),  rs3495  (G/  A) 

NP 

RAN  1-1 

rsll26680A/G 

NP 

rs55781031  (A/  G) 

NP 

rsl803274 (A/  G),  rs3495  (G/  A) 

NP 

RAN  1-2 

NP 

NP 

NP 

NP 

rsl803274 (A/  G),  rs3495  (G/  A) 

NP 

RAN  1-4 

rsll26680A/G 

NP 

rs55781031  (A/  G) 

NP 

rsl803274 (A/  G),  rs3495  (G/  A) 

NP 

RAN  1-5 

NP 

NP 

NP 

NP 

rsl803274 (A/  G),  rs3495  (G/  A) 

NP 

RAN  1-6 

NP 

NP 

NP 

NP 

NP 

NP 

RAN  1-7 

NP 

NP 

NP 

NP 

NP 

NP 

RAN  1-8 

NP 

NP 

NP 

NP 

NP 

rs77307863  (G/  T) 

RAN  1-9 

rsll26680A/G 

NP 

rs55781031  (A/  G) 

NP 

rsl803274 (A/  G),  rs3495  (G/  A) 

NP 

RAN  1-10 

NP 

novel  3  C/  T 

NP 

NP 

NP 

NP 

RAN  1-11 

NP 

NP 

NP 

NP 

rsl803274 (A/  G),  rs3495  (G/  A) 

NP 

RAN  1-12 

NP 

rsl799807A/G 

NP 

NP 

rsl803274 (A/  G),  rs3495  (G/  A) 

NP 

RAN  1-13 

NP 

NP 

NP 

NP 

NP 

rs77307863  (G/  T) 

RAN  1-14 

NP 

NP 

NP 

NP 

rs3495(G/A) 

rs77307863(G/T) 

RAN  1-15 

NP 

NP 

NP 

NP 

rsl803274 (A/  G),  rs3495  (G/  A) 

NP 

RAN  1-16 

NP 

NP 

NP 

NP 

rsl803274 (A/  G),  rs3495  (G/  A) 

rs77307863  (G/  T) 

RAN  1-17 

NP 

novel  1 C/  T 

NP 

NP 

rs3495(G/A) 

NP 

RAN  1-18 

NP 

rsl6849700  C/  G 

NP 

NP 

rsl803274 (A/  G),  rs3495  (G/  A) 

rs55849903  (C/  T) 

RAN  1-19 

NP 

NP 

NP 

NP 

NP 

rs77307863  (G/  T) 

RAN  1-20 

NP 

NP 

NP 

NP 

rs3495(G/A) 

NP 

RAN  1-21 

NP 

rsl799807A/G 

NP 

NP 

rsl803274 (A/  G)  rs3495  (G/  A) 

NP 
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Sample 

Exon  1 

Exon  2 

intron  2  (near  3'  ) 

Exon  3 

Exon  4 

3UTR 

RAN  1-23 

NP 

NP 

NP 

NP 

rsl803274  (A/  G) 

NP 

SOTG8-55 

NP 

NP 

NP 

NP 

NP 

NP 

SOTG8-56 

NP 

NP 

NP 

NP 

NP 

NP 

SOTG8-57 

NP 

novel  2  A/  G 

NP 

NP 

rs3495(G/A) 

NP 

SOTG8-58 

NP 

NP 

NP 

NP 

rs3495(G/A) 

NP 

SOTG8-59 

NP 

NP 

NP 

NP 

rs3495(G/A) 

NP 

SOTG8-60 

rsll26680A/G 

NP 

NP 

NP 

rsl803274 (A/  G);  rs3495  (G/  A) 

NP 

SOTG8-61 

NP 

NP 

NP 

NP 

NP 

NP 

SOTG8-63 

NP 

NP 

NP 

NP 

NP 

NP 

SOTG8-65 

NP 

NP 

NP 

NP 

NP 

rs77307863  (G/  T) 

SOTG8-66 

NP 

NP 

NP 

NP 

NP 

NP 

SOTG8-67 

NP 

NP 

NP 

NP 

NP 

NP 

SOTG8-68 

NP 

NP 

NP 

NP 

rsl803274(A/A) 

NP 

SOTG8-69 

NP 

NP 

NP 

NP 

NP 

NP 

SOTG8-70 

NP 

NP 

NP 

NP 

NP 

NP 

SOTG8-72 

NP 

NP 

NP 

NP 

NP 

NP 

SOTG8-73 

NP 

NP 

NP 

NP 

NP 

NP 

SOTG8-75 

NP 

NP 

NP 

NP 

NP 

NP 

SOTG8-76 

rsll26680A/G 

NP 

rs55781031  (A/  G) 

NP 

rsl803274(A/A) 

NP 

SOTG8-62 

NP 

NP 

NP 

NP 

NP 

NP 

SOTG8-64 

NP 

NP 

NP 

NP 

rsl803274 (A/  G),  rs3495  (G/  A) 

rs77307863  (G/  T) 

SOTG8-77 

rsll26680A/G 

NP 

rs55781031  (A/  G) 

NP 

rsl803274  (A/  G) 

NP 

N  P:  N  on  Polymorphic 
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Allele  frequencies  for  SN  Ps  were  calculated  when  appropriate  and  this  was  compared  to  the 
literature,  Table  4.  Overall,  this  small  ADF  sample  demonstrated  genotypefrequencies  similar 
to  AFD_EUR_PANEL  and  FlapMap-CEU,  both  Caucasian  panels  (Appendix  C).  For  detailed 
definitions  of  the  populations  shown  refer: 

http://ccr.coriell.orq/  Sections/  Collections/  NFIGRI/  hapmap.aspx?Pqld=266&coll=HG  and  a 

general  description  can  be  found  in  Appendix  D. 


Table4:  Allele  frequencies 


SN  P 

Population  (NCBI) 

ADF 

n=51 

Caucasian 
AFD  EUR  P 

ANEL 

H  apM  ap- 
CEU 

(European) 

H  aoM  ap- 

CEU 

African 
American 
AFD  AFR  P 

ANEL 

HapMap- 

YRI 

HaoMap- 

YR1 

Asian 

AFD  CHN 

PANEL 

HapMap- 

HCB 

HaoMap- 

HCB 

H  apM  ap- 
JPT 

H  apM  ap- 

1EI 

rsll26680 

A  0.058 

G  0.9411 

A:  0.083 

G:  0.917 

A:  0.058 

G:  0.942 

1.00 

1.00 

1.00 

1.00 

1.00 

rsl799807 

A. 0.980 

G. 0.200 

A  0.979 

G  0.021 

0.992 

0.008 

0.978 

0.022 

1.00 

1.00 

1.00 

1.00 

rsl6849700 

C  0.98 

G  0.02 

C  1.00 

1.00 

C  0.935 

G  0.065 

C  0.925 

G  0.075 

1.00 

1.00 

1.00 

rs55781031 

A  0.951 

G  0.049 

A  0.917 

G  0.083 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

rsl803274 

A  0.235 

G  0.765 

A  0.167 

G  0.833 

0.175 

0.825 

0.152 

0.848 

0.167 

0.833 

0.083 

0.917 

0.200 

0.800 

0.205 

0.795 

rs3495 

A  0.235 

G  0.765 

A  0.771 

G  0.229 

N/A 

A  0.341 

G  0.659 

N/A 

A  0.771 

G  0.229 

N/A 

N/A 

rs77307863 

A  0.09 

G  0.901 

A  0.056 

G  0.944 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/  A-  unknown,  no  frequency  data 


4.  Discussion 

BCFIE  iswidely  distributed  throughout  the  body.  TheBCHE  in  plasma  is  synthesized  in  the 
liver.  It  isa  biological  scavenger  against  organophosphorus  and  carbamatecompoundsand  has 
an  i  important  rol e  i  n  the hydrol ysi s  of  the  musd e  rel axant  sued  nyl chol  i  ne.  M  utati ons  in  BCHE 
gene  affect  the functi  on  of  thi  s  enzyme  I  ead i  ng  to  prol  onged  neu romuscul  ar  bl  ockade. 

This  study  successfully  screened  the 4  coding  regions  (and  surrounding  intronic  regions)  of 
BCHE  for  both  novel  and  known  polymorphisms  of  interest.  Overall,  35/  51individualscarry  a 
BCHE  polymorphism,  although  thefunctional  relevance  for  many  of  these  SN  Ps  is  unknown. 
FIRM  and  sequencing  procedures  revealed  12  different  genetic  polymorphisms  and  this 
included  the  K  variant  and  the  Atypical  gene,  along  with  4  novel  mutations. 

Besides  the  normal  phenotype,  the  most  frequent  variations  of  BCHE  are  the  Atypical  and  K 
variants.  Both  variants  have  been  associated  with  prolongation  of  muscle  relaxation  after 
administration  of  sued  nyl  choline  and  mivarcurium  [13].  In  Caucasian  populations  the  allele 
frequency  of  theAtypical  variant  is  about  2%,  whereas  the  allele  frequency  oftheK  variantis 
betw  een  16. 7- 17. 5%  w  i  th  a  h  omozy  go  u  s  i  n  ci  d  en  ce  of  1  i  n  ap  p  roxi  matel  y  63  i  n  d  i  v  i  d  u  al  s  [  13- 15] . 
The  K  variant  has  been  reported  to  be  the  most  common  variantin  any  population  and  in  vitro 
enzyme  activity  studies  have  demonstrated  a  33%  reduction  in  enzyme  activity.  In  theADF 
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population  examined,  both  the  Atypical  and  K  variant  were  identified.  The  K  variant  was 
observed  in  23individuals(21heterozygousand  2  homozygous),  representing  a  frequency  of  25 
of  102  alleles  (24%  allele  frequency).  This  allele  frequency  is  slightly  higher  than  that  found  in 
the  H  apM  ap  samples  of  randomly  selected  healthy  individuals  of  N  orthern  European  (17.5%) 
and  Asian  (20%)  decent.  At  present,  the  BCHE  activity  in  the  ADF  samples  has  not  been 
measured,  although  it  has  been  reported  that  both  heterozygous  and  homozygous  individuals 
with  theK  variant  have  reduced  serum  BCHE  activity  [10].  As  expected,  the  Atypical  variant 
was  found  at  a  much  lower  frequency,  with  only  2(1.8%)  heterozygous  individuals  carrying  the 
allele  which  is  consistent  with  frequencies  reported  for  Caucasian  populations. 

Previousstudieshaveshown  linkagedisequilibriumbetween  Atypical  and  K  variants  in  which 
the  K  vari  ant  w  as  fou  nd  i  n  89-96%  of  A  ty  p  i  cal  vari  ant  carri  ers  [9] .  Thi  s  has  show  n  to  be  the  case 
in  both  European  and  American  populations  however,  theAtypical  allele  has  been  shown  to  be 
rare  in  Japanese  [9, 10].  Genetic  investigations  have  also  provided  evidence  that  whilethe  K 
variant  alone  causes  moderate  prolongation  of  the  action  of  mivacuirum,  the  effect  is  more 
pronounced  in  patients  carrying  both  theAtypical  and  K  variants  [16].  Our  results  (although 
from  a  small  population)  support  this  as  both  samples  RANI- 12  and  RANI-21  in  which  the 
Atypical  variant  was  detected,  were  also  K  vari  ant  carri  ers. 

H  RM  analysis  identified  a  total  of  4  novel  BCHE  mutations,  3  in  exon  2  and  oneinthe3'UTRof 
BCHE.  Thesewereconfirmed  by  DNA  sequencing.  The  mutations  named  "novel  2"  and  "novel 
3"  describe  the  nucleotide  substitution  that  leads  to  an  amino  acid  change  lle384Val  and 
Ala212Val,  respectively.  Further  studieswill  comparethe  protein  level  sand  activity  of  BCHE  in 
theseindividualstodetermineany  possible  functional  consequence  of  these  polymorphisms. 
Thethird  mutation  "novel  3"  wasasynonymousbasepairchange,aC-T  at  position  1244and  is 
uni  i  kely  to  have  any  effect  on  the  expression  or  activity  of  6 CH  E .  Thefi  nal  mutati  on  found  i  n 
BCHE  is  located  in  the  3'UTR,  a  C-A  substitution  at  position  2366  of  the  cDN  A,  may  affect 
mRNA  stability.  Futureworkwill  comparetheB  CHE  protein  levelsofthisindividual  to  control 
samples. 

I  n  cond  usi  on,  the  outcomes  of  thi  s  study  are: 

•  FIRM  is  ableto  effectively  screen  for  and  identify  polymorphisms  in  theBCHE  gene. 

•  This  approach  detected  eight  known  and  four  novel  polymorphisms  in  a  small  ADF 
population. 

•  Two  of  the  most  commonly  studied  functional  relevant  polymorphisms  theK  variant 
and  A  ty  pi  cal  vari  ant  were  i  d  enti  f  i  ed  i  n  thi  s  su  bset  of  i  nd  i  vi  d  ual  s  f  rom  the  A  D  F . 

•  Comparison  to  SNP  databases,  revealed  the  genotype  frequencies  of  the  8 
polymorphisms  identified  were  similar  to  those  found  in  Caucasian  populations. 

•  Linkagestudiesin  European  and  American  populations  identified  theAtypical  alleleis 
often  carried  with  the  K  variant.  Our  data  support  these  findings. 
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These  findings  are  significant  as  they  provide  a  cheap  and  rapid  method  of  obtaining  prior 
knowledge  of  possible  hypersensitivity  to  sued nyl choline  and  various  organophosphorus 
compounds  in  members  of  the  AD  F. 

Further  work  will  involve  the  generation  of  a  larger  data  set  to  allow  a  more  accurate 
determi  nati  on  of  al  I  el  efrequenci  es  of  polymorphisms  of  i  nterest  to  the  A  DF.  The  measurement 
of  BCHE  protein  and  activity  in  serum  samples  will  allow  the  determi  nati  on  of  the  functional 
rel  evance  of  the  novel  vari  ati  on  found . 
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Appendix  A:  BCHE  polymorphisms 


It  should  be  noted  that  the  naming  convention  of  SN  P  location  has  changed  for  bche  over  theyears.  Initial  naming  of  SN  Ps  used  their  position 
rel  ati  veto  N  -terminus  of  the  mature  protein,  listed  as  the  alternative  name  [17].  Current  existing  convention  isto  number  the  effected  ami  no  acid 
from  the  signal  sequence  rather  than  the  N-terminus  of  the  mature  protein  (effect  column).  Where  available,  database  reference  numbers  (rs 
numbers) areprovided.  Nudeotidechangeisrelativetothefull  length cDNA,  referenceNM_000055.2orthegenomicsequenceNT_005612.16if 
theSNP  is  located  within  an  intron. 


TableAl.  Description  cfthemain  polymorphism  identified  in  thelita^atureand  pubic  databases.  Wha'e  avail  able  theSNP  ID,  al  ta^nati  ve  name 
nucleotide  change  and  dfect  (if  any)  on  the  protan  islistad. 


dbSNPID 

Alternative  name 

Nudeotide  change 

Effect 

Functional  studies 

Reference 

In  vivo 

In  vitro 

rs55950599 

N/A 

A176T 

T6S 

Unknown 

Unknown 

dbSNP  build  132 

Silent  2 

BCH  E*FS6 

A  259- 

Frame  shift  at 
codon  6 

Silent  phenotype 

deletion  of  A  at  position 
259  of  full  length  cDNA 

N  o  protei  n  expressed 

[11] 

rsl 14063013 

N/A 

A227G 

123V 

Unknown 

Unknown 

dbSNP  build  132 

rsll6047990 

N/A 

A279G 

K40R 

Unknown 

Unknown 

dbSNP  build  132 

rs56309653 

BCHE*24M 

C315T 

T52M 

Unknown 

Reduced  activity  of 
enzyme 

[10] 

rsll4355070 

N/A 

C324A 

A55D 

Unknown 

Unknown 

dbSNP  build  132 

BCHE  281 

T326A 

F56I 

Suggested  pathogenic 
variant 

BChE  deficiency 

[12] 

rsll6097205 

BCH  E*33C 

A342G 

Y61C 

Silent  phenotype 

Normal  protein  levels  in 
plasma  enzymatically 
inactive  in  plasma. 

[11] 

BCH  E*37S 

C353T 

P65S 

Silent  phenotype 

Protein  expressed  at  low 
levels  in  culture.  Poor 
activity  1-2% to  BZ,  BTC, 
ACT,  PTA 

[11] 

rs75995351 

N/A 

C373A 

F71L 

Unknown 

Unknown 

dbSNP  build  132 

UNCLASSIFIED 


15 


DSTO-TR-2503 


UNCLASSIFIED 


dbSNPID 

Alternative  name 

N  ud  eoti  de  change 

Effect 

Functional  studies 

Reference 

In  vivo 

In  vitro 

rsl799807 

Atypical  BCHE 
BCHE*70G 

A453G 

D98G 

Silent 

Sued  nylcholine  sensitivity 
Pyridostigminesensitivity 

Dibucaine  resistant 
Predicted  to  have 
increased  sensitivity  to 
positively  charged 
organophosphates,  no 
diff  with  uncharged,  eg. 
paraoxon  and  DFP 

[18] 

BC  H  E*75R 

G467C 

G103R 

reduced 

Unknown 

45%  of  activity  BTC  as 
compared  to  wt  in  293T 
kidney  cells 

[19] 

BC  H  E*90D 

A514C 

E118D 

Silent  phenotype 

Protein  rapidly 
degraded,  no  activity  in 
293T  cells 

[19] 

BCHE*N96f 

A53CTT 

N 124Y 

Unknown 

Unknown 

Australian  specific 
[12] 

rs3732880 

BCH  E*100S 

C542T 

P128S 

Silent  phenotype 

Reduced  protein  activity 

[10] 

BCH  E*106fs 

A  insertion  at  562 
of  cDNA 

N134FS 

Silent  phenotype 

Insertion  of  a  stop 
codon,  23  ami  no  acids 
downstream 

Australian  specific 
[12] 

Silent  7 

BCH  E*115D 

G591A 

G115D 

Silent  phenotype 

[20] 

rs55703238 

si  lent-1  BCH  E*FS117 

-594T 

G117fs 

Silent  phenotype 

Deletion  of  aT  causes  a 
shift  in  thereading 
frame  from  Glyll7 
where  GGT(Gly)- 
GGAG  (Gly  and  one 
base)  to  a  new  stop 
codon  at  position  129 

N  o  protei  n  expressed 

[8] 

BCH  E*125F 

A619T 

L153F 

Silent  phenotype 

Protein  expressed  at  low 
levels  in  culture.  Poor 
activity  1-2% to  BZ,  BTC, 
and  PTA.ACT 
improved 

[11] 
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dbSNPID 

Alternative  name 

Nudeotide  change 

Effect 

Functional  studies 

Reference 

In  vivo 

In  vitro 

H  variant 

BCH  E*142M 

G669A; 

V170M 

Sued  nylchol  i  ne  sensitivity 

A  p  proxi  matel  y  10% 
reduced  concentration 

[19] 

BC  H  E*170E 

T745G 

D198E 

Silent  phenotype 

Protein  expressed  at  low 
levels  in  culture.  Poor 
activity  1-2% to  BZ,  BTC, 
ACT,  PTA 

[11] 

rsll4706984 

BCH  E*184V 

SC  variant 

C795T 

A  212V 

Unknown 

Plasma  sample  showed 
decreased  ability  to 
hydrolyse 
sued  nyl  choline 

[21] 

BCH  E*198G 

A836G 

S226G 

Silent  phenotype 

Normal  protein  levels, 
enzymatically  inactivein 
plasma, 

[11] 

BCH  E*199V 

C841T 

A 199V 

Unknown 

Same  protein  levels  and 
enzyme  activity  as  wt  to 
BTC  in  293T  kidney  cells 

[19] 

BCH  E*201T 

G847A 

A229T 

Silent  phenotype 

Normal  protein  levels, 
enzymatically  inactivein 
plasma 

[11] 

rs28933389 

fluoride- 1 
BCHE*243M 

C972T 

T271M 

Partial  sued nyl choline 
sensitivity 

Resistant  to  inhibition  by 
0.050  mM  sodium 
fluoride  in  the  in  vitro 
assay 

[7] 

rsl6849700 

G1009C 

R283D 

Unknown 

Unknown 

dbSNP  build  132 

rsll5624085 

A1044G 

K295R 

Unknown 

Unknown 

dbSNP  build  132 

rsll4991020 

A1094T 

T312S 

Unknown 

Unknown 

dbSNP  build  132 

rsll6607681 

C1119T 

P320L 

Unknown 

Unknown 

dbSNP  build  132 

rs!04893684 

BCHE*335P 

T1164C; 

L335P 

silent 

Silent  phenotype 

No  protein  in  plasma 

[22] 

BCH  E*271X 

G1055T 

E299X 

Silent  phenotype 

A  bsence  of  protei  n  and 
activity  in  human 
plasma 

[11] 
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dbSNPID 

Alternative  name 

Nudeotide  change 

Effect 

Functional  studies 

Reference 

In  vivo 

In  vitro 

BCH  E*328D 

C1227A 

A355D 

Predicted  silent  phenotype 

Unknown 

[23] 

BCHE*K355insALU 

1306- 1320  Alu 
insertion 

K383AIU 

Silent  phenotype 

Insertion  of  an  ALU 
element  at  codon  383 

[24] 

BCHE*FS315 

ACC-AACC  at 
position  1187 

T343fs 

Silent  phenotype 

A  insertion  at  position 
1187  of  cD  N  A . 

Frameshift  at  codon  343 

[25] 

rsll5129687 

BCH  E*365R 

G1337C 

G393R 

Silent  Phenotype 

no  protein 

[25] 

rs35979453 

T1357- 

F399X 

Unknown 

frameshift 

dbSNP  build  132 

BCH  E*R386C 

C1400T 

R413C 

Unknown 

Unknown 

[12] 

Australian  specific 

rs28933390 

fluoride-2 

F-2;  BCHE*390V 

G1413T 

G418V 

Patients  have 
sued  nylchol  i  ne  sensitivity 

Fluoride  resistant 

[26] 

BCH  E*418S 

T1497V 

F446S 

Silent  phenotype 

[27] 

BCHE424X 

C1515T 

R452X 

Silent  phenotype 

Studies  suggest  non 
functional  protein 

Australian  specific 
[12] 

BCH  E*E460K 

G1622A 

E488K 

Unknown 

Unknown 

[12] 

BCH  E*465X 

A1637T 

R493X 

Silent  phenotype 

Premature  stop  codon 

593 

[27] 

rsll5017300 

C1652T 

R498W 

Unknown 

Unknown 

dbSNP  build  132 

BCH  E*471R 

T1655 

W499R 

Silent  phenotype 

Protein  expressed  at  low 
levels  in  culture.  Poor 
activity  1-2% to  BZ,  BTC, 
ACT,  PTA 

[11] 

J  variant 
BCHEE497V 

A 1734T ;  exon  3 

E524V 

Sued  nylchol  i  ne  sensitivity 

33%  reduced 
concentration 

[28] 

BCH  E*500X 

T1745A 

528X 

Silent  phenotype 

A  bsence  of  protei  n  and 
activity  in  human 
plasma 

[11] 

rs77586249 

T1773C 

I538T 

Unknown 

Unknown 

dbSNP  build  132 

BCH  E*515C 

C1787T 

R543C 

Silent  phenotype 

[27] 
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dbSNPID 

Alternative  name 

Nudeotide  change 

Effect 

Functional  studies 

Reference 

In  vivo 

In  vitro 

BC  H  E*518L 

A1797T 

Q546L 

Silent  phenotype 

Protein  expressed  at  low 
levels  in  culture.  Poor 
activity  (1-2%)  to  BZ, 

BTC,  ACT,  PTA 

[11] 

rsll6333990 

C1815T 

S552L 

Unknown 

Unknown 

dbSNP  build  132 

rsll4166903 

T1851C 

I564T 

Unknown 

Unknown 

dbSNP  build  132 

rs56325145 

BCHE*I3E4-14C 

T685-14C 

165491308 

In  linkagewith  BCHE*115D 
(known  silent  mutation) 

Unknown,  nucleotide 
substitution  in  intron  3, 

14  nucleotides  upstream 
of  exon  4 

[23] 

rsl803274 

K  variant 

BCH  E*539K 

G1859A 

K567T 

Succinylcholinesensitivity 

66%  reduced 
concentration 

[9] 

rs2668207 

165493724* 

C/T 

Association  with  lower 
chol  i  nesterase  1  evel  s 
independent  of  pesticide 
exposure 

Located  within  intron  2, 
43-9.5%  effect  on 
activity 

[29] 

rs2048493 

165544302 

C/G 

Association  with  lower 
chol  i  nesterase  1  evel  s 
independent  of  pesticide 
exposure 

Located  within  intron  2, 
43-9.5%  effect  on 
activity 

[29] 

DN :  the  percent  inhibition  of  activity  caused  by  0.03umol.L  dibucaineBZ:  benzoylcholine;  BTC;  butyrylthiocholine,  ACT;  acetylthiocholine,  PTA;  phenylthioacetate 
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Appendix  B:  Sequencing  results  of  HRM  Analysis 

This  appendix  presents  an  example  of  the  sequencing  data  obtained  for  each  unique  SNPs 
identified  in  this  study.  In  most  cases  a  SN  P  was  deemed  confirmed  if  a  second  independent 
PCR  was  performed  and  the  SN  P  was  observed  both  the  sense  and  anti-sense  directions,  or 
using  2  different  primers. 

Exon  1  fragment  1  rsll26680  (G  /A ). 

SOTG8-76(G/  A) 


Exon  2  Fragment  3  rl799807  (G/A) 

RAN  1-21  (G/A) 
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i  f  j  x  r  x  r  r  iiwiiiji  r  r  r  ■  r  i  1  1 1  j  §  1  *  1  iji  1  i  ■■ij  biiijiIbwii  x  x  1  i  fiii  x  x  r4T  f  f  i  x  ■  t  r  r  r  ■  j  ■  if  i 

wA--/Vw'/y /  Y  A  /yW'-  .vw^v/ . vv^VwA  / VYvvvvYvw'A 

1 

d 

r  11  f  X  X  1  1  t  f  f  f  1  f  1  f  1  f  f  f  1  1  iiliil]ll  1  1  f  f  1  f  1  f  1  1  X  X  f  f  I  X  1  1  X  Xf  XfjiX  X  f  J  fii  1  I  1  f  f  f  I 

Y\mAa/WaAAAaaAaA/W\aaaaaaaM/\AaaAaAM 

*■- 

i 

d 

if  (  r  r  t  r  t  1 1 1 

uAaAaAaAJ 

.  -  . ..  .  _ 
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i- 
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- 1 - 1 - - - — — — - A - - — * — - -  j\\d  W  W  iA/\A  J 

Exon  2  Fragment  4  N  ovel  (C/T) 


RAN  1-10  (C/  T) 


J 


Hi  e  f  i  tit  i  n  i  fin  ttiH  nmmim  tn  t?  h  i  i  t  h  tit  mt  h  tan  1 1  nifmii  u  t  a 
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a. 


— - 


AvWWyyvWW 
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S  6  B  r  T  C 


Exon  2  Fragment  6,  rsl6849700  (G  /C ) 
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RAN  1 18 (G/  C) 


Aiwk 


'Ac 
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- - 
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lAAAMAAAAMA^AAAAAMAAMAAAAAA/^Al/\AA/\A^/^l\AAA|\/^ 


a/WvLaaaAa 


•  /\  A  . .  -'"•-/wy  'i  w  .  A.  A. . .  /\ 


Vi/VA~  A.  Ail 


Exon  2  Fragment  7.2  N  ovel  1  (C/T) 

RAN  1-17  (C/T)  ^ 


tiN-»  Nfi  HHiti. 


■H-  it  f  f  -H  Hrt-t-t  -t-H-H  M  i  H-  f  t~H-t  i 

Aa/AAA/'A-'/V-.--  JV:-- -•■A--.  A"  ■  -Aaa  - 


M 


^\Aa/AaaAVV^\aaA/w,^  \ 


A 


f'H-HfjNifHH-HH-l  f  f-H  tti  t-tt-i-H-i  i  i-f  Hfj  HHirHfH  »  <-i-ii-  i-hf- 1  +-(■■«■+  il  H-it-f  i 


iM/VkAil 


'AJAMNWvMm 


qAAMz 


A 


Exon  2  Fragment  7.2  N  ovel  2.  (A/G ) 

SOTG8-57  (A/  G) 


1 


n  H  h  i  f  f  i  i  h  ( t  f  m  f  n  ^  t  m  u  m  n  i  >  f  f  <  j  u  j  1 1  n  m  >  m  n  m  n  <  n  i  n  >  >  n  n  aj  i  i 

.aA/v'AaaaaaAaAQ.AAaMAaa/1A&aAA/Ww\MAAAa1AA 


WWlMAAAAAAAMAAMMAAMAA/lAAi 


Exon  3  Fragment  11  rs55781031  (A/G ) 
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Exon  4  Fragment  13.1,  rs!803274  (G/A).  All  Allele  combinations  detected. 


RAN  1-13  (G/G). 


1 

: 

WM'iMNtil 

W/wkVvWw A/WWV\ 

■  ' 

i 

/WW\MA 

V  i/VvWV 

1  1 

il;v.\A l|-1'  '■ 

SOTG8-47  (G/  A) 


SOTG8-76(A/  A) 


Exon  4  Fragment  14  rs3495  (G/  A ) 


SOTG8-59(G/  A) 


.1  Hi  J  H  iii-Mi  UH-t  H-t+i  mtNHHMn  TWt  <1  ii  it-  tii  +  t 


ft 


-ti  f  fi t-tt-t-i 


7~V~7 - 7“7 - 7  7 — 7" 


31 


± 


_ 


Mi/WWV 
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Exon  4fragment  15,  rs77307863  (G/  T).  All  all  Recombinations  found 


RAN  1-2  (G/G) 


RAN  1-2- F3 oontrol 0 F05  Fragment  base  #268.  Base  212  o*  238 

G  GTC  TC  T  TTGTTTTTTRTRTRRRGGCRRGP 

)AaaA/\MaAAaAAAA/V\a^aaAaaa 

,aaaA 

1TCT  TRRRRT  GR  GRGG  C  TCflGCN 

RAN  1-16  (G/  T) 

j _ 

from  RANI  16  r*77307863  (G7T) 

t  r  c  crrrrr  r  r  e  r  r  r  i  r  i  ff  r  a  r  A  a  c  c  P  fv  ff  e 

MAAAAAMAMAAAAAAfWWuA/WAAl 

Tirol 

a/WV 

1  Cl  1HHHRI  bRuRuG  G  I  C  H  b  L  N 

aAaaAaaAaMaAaA/nAMjAaIIJ 

SOTG8-52  (T/  T) 

BGIC1C1  1161  T1I11AIA1AAA66CA  a"*g’"/ 

l»  lr»}f-*rt  n>t»  r,m 
A  A  B  A  ; 

Him  «1  fit  241 

1  C  1  IAAAA1  6  A  G  A  G  (i  C  1  CAli  C 

GGirinrilf.1  IIIlIHtRIRRHubCRRGH 

A  aaAaAA  A  AAaAAAAaAaAM  a  aAamAA/ 

H  H  1  R  f 

A  / 1 A  A  f 

rffiHRRBiGHCHCc  r  r  f  h  c  f  >3 

aAaaAaaAaAAaAaA/vAM^aA 

Exon  2  Fragment  15  rs55849903  (C/  T) 


RAN  1-18  (C/  T) 


4-4  *  -  ■.  J«I  1.  -• 

-S'  J  -*3 

■«i**  HI. C  'i»  * r -*. •  -  r  *  mh  - 

7  4  4  a  f  4 

\A^vvv\^aa/\AA^/v\aMAaMV\^^  a/  4aV.  .. 

1 

C  f  i  f  i  6  f  1  t  1  1  It  Cl  C  U  (  C  1C  1  h  1  1  II  Chi  I  k  ({ h  c  rc«-  lalflllkNCCfhl  A  (  1  1  1  (  1  ft  ft  »  i  (III  I  ft  O'  ft  f  ft  f  1  C  ft  ft  ■' 

AWWWvv-  <  vw\ /vvA- , / mm  ,\a  ,  vwww  m/W  .VWt^ i/Vv,  aAA/iv:;;.  .7,  'AaaaaA 

Exon  2  Fragment  15  novel  4  (C /  T) 
SOTG8-46(C/  T)  ^ 


I0>»  •  -j  U.wf  >IHM  mm«'¥  ln<  Mt  .1  iV 

NNHimuHNUfiinf 


■rut  t-H  it  t  tut  ft  ut  i  &  u  i  tt  i  iWriWrrrfH 

An "  ‘ A  ^/WWV\MMAA/W\AAAA/\A<vMA/MMfflA/ 
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Appendix  C:  Genotype  frequencies 

This  appendix  describes  the  genotypefrequencies  of  theSN  Ps  identified  in  this  study  as  compared  to  external  databases. 


Tabled:  Thegaiotypefrequeidts  cf  SNPs  cf  interested  as  compared  to  theNCBI  and  db6N  P  132  database 


SNP 

Population  (NCBI) 

ADF 

n=51 

Caucasian 

AFD  EUR  PANEL 

HapMap-CEU 

(European) 

HapMap-CEU 

African  American 

AFD  AFR  PANEL 

HapMap-YRI 

HapMap-YRI 

Asian 

AFD  CHN  PANEL 

HapMap-HCB 

HapMap-HCB 

HapMap-JPT 

HapMap-IPT 

rsll26680 

0.167  (A/ G) 

0.833  (G/G) 

0.017  (A/ A) 

0.050  (G/ A) 

0.900  (G/G) 

1.000  (G/G) 

1.000  (G/G) 

1.000(G/  G) 

1.000  (G/G) 

1.000  (G/G) 

rsl799807 

0.958  (A/ A) 

0.042  (A/ G) 

0.983  (A/ A) 

0.017  (A/ G) 

1.000  (A/ A) 

1.000  (A/ A) 

1.000  (A/ A) 

1.000  (A/ A) 

1.000  (A/ A) 

rsl6849700 

1.000  (C/C) 

1.000  (C/C) 

0.870  (C/C) 

0.130  (G/  C) 

0.850  (C/C) 

0.150  (G/  C) 

1.000  (C/C) 

1.000  (C/C) 

1.000  (C/C) 

rs55781031 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

rsl803274 

0.042  (A/ A) 

0.250  (G/ A) 

0.708  (G/G) 

0.017  (A/ A) 

0.317  (G/ A) 

0.667  (G/G) 

0.304  (A/ G) 

0.696  (G/G) 

0.317  (A/G) 

0.683  (G/G) 

0.167  (A/G) 

0.833  (G/G) 

0.022  (A/ A) 

0.356  (A/G) 

0.622  (G/G) 

0.067  (A/ A) 

0.267  (A/G) 

0.667  (G/G) 

rs3495 

0.583  (A/ A) 

0.375  (A/  G) 

0.042  (G/G) 

0.045  (A/ A) 

0.591  (A/G) 

0.364  (G/G) 

0.542  (A/ A) 

0.458  (A/G) 

0.000  (G/G) 

rs55849903 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

rs77307863 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 
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Appendix  D :  Population  descriptions 


Table  Dl:  Summary  of  populations  and  number  of  individuals 


Population 

Approved  individual  DNA  samples 

Approved  Individual  cell 
cultures 

Yoruba  in  Ibadan,  Nigeria  [NRI] 

220 

220 

Han  Chinesein  Beijing,  China 
[CHB] 

162 

162 

JapaneseTokyo,  Japan,  [YPT] 

131 

131 

CEPH  collection  [CEU] 

180 

180 

A  brief  description  of  the  Hap  Map  (HapMap)  populations  is  provided  below.  Additional 
information  can  be  found  at  http:/  /  hapmap.ncbi.nlm.nih.gov/  citi  nqhapmap.html  ,en 

Yoruba  in  Ibadan,  Nigeria  (YRI) 

These  samples  were  collected  in  a  particular  community  in  Ibadan,  Nigeria,  from  individuals 
who  identified  themselves  as  havi  ng  four  Yoruba  grandparents.  I  nd  udi  ng  the  name  of  the  city 
and  cou  ntry  w  here  these  parti  cu  I  ar  Y  oruba  sampl  es  were  col  I  ected  al  so  rei  nforces  the  poi  nt  that 
the  sample  set  does  not  necessarily  represent  all  Yoruba  people,  whose  population  history  is 
complex.  These  samples  should  not  be  described  merely  as  "African,"  "Sub-Saharan  African," 
"West  African,"  or  "N  igerian,"  si  nee  each  of  those  designators  encompasses  many  populations 
with  many  different  ancestral  geographies. 

Japanese  in  Tokyo,  Japan  (JPT) 

These  sampl  es  were  col  I  ected  i  n  theT  okyo  metropol  i  tan  area,  from  peopl  e  who  came  from  (or 
whose  ancestors  came  from)  many  different  parts  of  Japan.  These  samples  arerepresentativeof 
the  majority  population  injapan,  as  it  included  sampl  es  from  peopl  ewhosegrandparents  were 
all  from  Japan. 

Han  Chinese  in  Beijing,  China  (CHB) 

Samples  werecollected  from  individuals  I  iving  in  the  residential  community  at  Beijing  Normal 
University  who  were  self- identified  as  having  at  least  three  out  of  four  Han  Chinese 
grandparents.  Although  individuals  of  Beijing  University  were  from  many  different  parts  of 
China,  this  set  of  samples  was  not  drawn  to  be  representative  of  all  Han  Chinese  people. 

CEPH  (Utah  Residents  with  Northern  and  Western  European  Ancestry)  (CEU) 

These  samples  were  collected  from  people  living  in  Utah  with  ancestry  from  northern  and 
western  Europe.  Theterm  "CEPH  "standsfortheCentred'Etudedu  PolymorphismeHumain, 
the  organisation  that  collected  these  samples  in  1980.  Because  the  importance  of  precision  in 
assigning  group  membershipto  prospective  donors  based  on  ancestral  geography  was  not  well 
appreciated  in  1980,  itisundear  how  accurately  these  samples  reflect  the  patterns  of  genetic 
variation  in  peoplewith  northern  and  western  European  ancestry.  These  samples  should  not  be 
described  as  "European,"  nor  seen  as  representing  people  with  ancestry  from  other  parts  of 
Europe  (eg.,  southern  or  eastern  Europe).  The  samples  also  should  not  be  described  as 
"C  au  casi  an , "  a  ter  m  that  car  ri  es  raci  al  overtones,  an  d  that  tech  n  i  cal  I  y  refers  on  I  y  to  peop  I  e  f  rom 
the  area  between  the  Black  and  Caspian  seas 
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The  "panel"  populations  are  described  in  detail  at: 

http://ccr.corieil.orq/  Sections/  Collections/  NHGRI/  hapmap.aspx?Pqld=266&coll=HG 

AFP  EUR  PANEL 

These  24  samp  I  es  from  the  Cori  el  I  Cel  I  Repository  are  pri  mari  I  y  of  European  A  meri  can  descent. 
Samples  consist  of  23  unrelated  CEPH  parents  selected  by  the  SeattleSN  Ps  Program  for 
GenomicApplications,  pi  us  one  sample(N  A 17201)  from  Cori  el  I  'shuman  variation  panel  of  50 
Caucasians  (HD50CAU). 

AFP  CHN  PANEL 

These  24  samples  of  Chinese  descent  from  the  Cori  el  I  Cell  Repository  were  selected  from  the 
hi  an  People  of  Los  Angeles  Panel  of  lOO(FIPlOOCFII). 

AFP  AFR  PANEL 

These  23  samples  of  African  American  descent  from  the  Cori  el  I  Cell  Repository  were  selected 
from  the  human  variation  panel  of  50  African  Americans  (FI  P50A A). 
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